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The Large Area Telescope (LAT), Fermi's main instrument, is providing a new view of the local energetic pulsar 
population. In addition to identifying a pulsar origin of a large fraction of the bright unidentified Galactic 
EGRET sources, the LAT results provide a great opportunity to study a sizable population of high-energy 
pulsars. Correlations of their physical properties, such as the trend of the luminosity versus the rotational 
energy loss rate, help identify global features of the gamma-ray pulsar population. Several lines of evidence, 
including the light curve and spectral features, suggest that gamma-ray emission from the brightest pulsars 
arises largely in the outer magnetosphere. 



1. INTRODUCTION 

Since the launch of the Fermi satellite and dur- 
ing its first six months of data taking in survey 
mode the pulsar population discovered in the 
gamma-ray sky is grown as far as 46 objects. All 
their features and a first study of the new popu- 
lation are summarized in The First Fermi Large 
Area Telescope Catalog of Gamma-ray Pulsars 
lAbdo. A. A., et al.. for the Fermi LAT collaboration! 
2009], Fermi data were analyzed to search for 
pulsations using 762 contemporaneous ephemerides 
obtained from radio telescopes and 5 from X-ray 
telescopes. A group of 218 of these ephemerides come 
from a list of good candidate gamma-ray pulsars 
selected with spin-down power greater than 10 34 
erg s _1 . The rest of the ephemerides is a sample 
of 544 pulsars from nearly the entire P — P plane, 
reducing in this way the possible bias of the LAT 
pulsar searches. In total the LAT detected 29 radio- 
selected puls ars, 5 of them are EGRET pulsars, and 
8 are MSPs [Abdo et al\ l2009d| . The L AT has also 
discov ered 16 gamma-selected pulsars [Abdo et~al\ 
2009a]. The discovery of these pulsars, as well 
as the determination of their timing ephemerides 
come directly from LAT gamma-ray data. Including 
Geminga, the new population counts 17 gamma- 
selected pulsars. A campaign to search for their radio 
counterpart started soon and three of them have 
alread y been found Camilo et al\ l2009j j [Abdo et~cd\ 



are true ra dio-quiet gamma-r ay pulsars, as in the case 
of CTA1 jAbdo et all l2008j and PSR 1836+5925, 
that were deeply searched for radio pulsations in the 
past. With this large population of pulsars statistical 
studies on their features can be performed, that 
allow to better understand their physics and emission 
mechanisms. 
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Figure 1: Phase difference A between the gamma-ray 
peaks, versus the phase lag 5 between the main radio 
peak and the nearest gamma-ray peak. Pulsars without a 
radio detection are plotted with 5 = 0. With present 
light curves we cannot generally measure A < 0.15; 
objects classified as single-peaked are plotted with A=0. 
Two such objects, both MSPs, are off the plot at S > 0.8. 
Blue squares: gamma-ray-selected pulsars. Red triangles: 
millisecond gamma-ray pulsars. Green circles: all other 
radio loud gamma-ray pulsars. 
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Figure 2: Value of the exponential cutoff Soutofi versus 
the magnetic field at the light cylinder, £?lc- The 
histogram of E cuto g values is projected along the 
right-hand axis. Blue squares: gamma-ray-selected 
pulsars. Red triangles: millisecond gamma-ray pulsars. 
Green circles: all other radio loud gamma-ray pulsars. 
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Figure 3: Photon index V versus the rotational energy 
loss rate, E. The histogram of the photon indices is 
projected along the right-hand axis. Blue squares: 
gamma-ray-selected pulsars. Red triangles: millisecond 
gamma-ray pulsars. Green circles: all other radio loud 
gamma-ray pulsars. 



2. TIMING AND SPECTRAL FEATURES 

Deep analysis were performed on all the pulsars 
detected by Fermi. Here we summarize and discuss 
their results. A detailed description of these analysis 
can been found in the first Fermi pulsar catalo 



lAbdo. A. A. . et al.. for the Fermi LAT collaboratio: 
2009]. 

Long base ephemerides that cover the first six 
months of LAT sky survey were evaluated for each 
one of the detected pulsars. For the radio loud pul- 
sars they are built joining the radio data of several 
radio telescopes. All the telescopes involved in this 
work are Parkes, Nancay, Westerbork Synthesis radio 
telescopes, Arecibo, Green Bank Telescope, and the 
76-m Lovell radio telescope at Jodrell Bank. For the 
gamma-selected pulsars the timing solutions are eval- 
ua ted using the LA T-data with a technique described 



Rav et qfJl2009j . The LAT observation is divided 



in several time windows of 10-20 days depending on 
the brightness of the source; the data in the windows 
are analyzed in order to define a pulse time of arrival 
(ToA) for each of them, that are finally used to build 
the time solution over the entire period of observation. 

The light curves obtained folding the photons 
with the calculated ephemeris were analyzed and 
parametrized. As first, they have been classified for 
their peak multiplicity, as single or double peaked. In 
the latter case the distance in phase of the two gamma- 
ray peaks (A) was measured. Most of the radio loud 
gamma-ray pulsars have a radio profile as simple as 
a single peak, or it is always possible in general to 
clearly define a main radio peak. The phase delay of 
the first gamma-ray peak respect to the main radio 
peak (S) was evaluated. 

Most of the pulsars show two dominant, relatively 



sharp peaks, suggesting the caustics from the edge of 
a hollow cone. When a single peak is seen, it tends 
to be broader, suggesting a tangential cut through an 
emission cone. This picture is realized in the Outer 
Gap and the high altitude portion of the Slot Gap 
models. As well, these models predict the correlation 
observed in figure [1] in the radio lag (S) versus the the 
gamma-ray peaks separation (A). 

The spectrum energy of all the LAT gamma-ray pul- 
sars is well fitted by a power-law with an exponential 
cut-off. 
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The significance of the exponential cut-off respect 
of a simply power-law model has been evaluated with 
the likelihood ratio test. For some pulsars the sig- 
nificance is below 5cr, and the cut-off energy (E c ) is 
poorly determined for those few with the significance 
lower than 3c. 

In figure is plotted the cut-off energy versus the 
magnetic field at the light cylinder (Blc)- 
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Along the right axis of figure [3] is shown the his- 
togram of the cut-off energies (E c ), that have a range 
of only about a decade, from 1 to 10 GeV, and are 
weakly correlated with the magnetic field at the light 
cylinder (B LC )- 

The fact that this features are common to all the 
different types of pulsars strongly implies that the 
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Figure 4: Galactic plane pulsar distribution (polar view). The star represents the Galactic center. The two circles 
centered at the Earth's position have radii of 3 kpc and 5 kpc. For pulsars with different possible distances, the nearer 
values are used. Blue squares: gamma-ray-selected pulsars. Red triangles: millisecond gamma-ray pulsars. Green 
circles: all other radio loud gamma-ray pulsars. Black dots: Pulsars for which gamma-ray pulsation searches were 
conducted using rotational ephemerides. Gray dots: Known pulsars which were not searched for pulsations. 



gamma-ray emission originates in similar locations in 
the magnetosphere relative to the light cylinder. Such 
a correlation of E c with B^c is actually expected in 
all outer magnetosphere models where the gamma- 
ray emission primarily comes from curvature radiation 
of electron s whose acceleration is balanced by radia- 
tion l o sses iMuslimoy and Harding|[2004 1 Zhang et all 
l2004j |Hirotanill2008j [12], 



In figure [3] is plotted the spectral index versus the 
rotational energy loss rate (E). As shown in the his- 
togram along the right axis, the spectral indices are 
distributed around ~ 1.5 and there is a general trend 
for the young pulsars to show a softer spectrum at 
large E. This may be indicative of higher pair multi- 
plicity, which would steepen the spectrum for the more 



energetic pulsars, either by steepening the spectrum of 
the curvature radiation-generating primary electrons 
[Romani 1996] or by inclusion of an additional soft 
spec tral component associated with robust pair form a- 
tion [Takata and Chanell2007t [Harding et 0/1120081 ]. 



3. THE LOCAL POPULATION AND THE 
PULSAR LUMINOSITY 

To evaluate the Galactic pulsar distribution as well 
as their luminosity, a good knowledge of the distances 
is needed. The annual trigonometric parallax method 
give the most precise estimate of the distance, but 
unfortunately it is available only for few nearby pul- 
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Figure 5: Gamma-ray luminosity L 7 versus the rotational energy loss rate E. Dashed line: L 7 equal to E. Dot-dashed 
line: L 7 proportional to the square root of E. L 7 is calculated using a beam correction factor fo = 1 for all pulsars. 
For the Crab is also plotted the total high energy luminosity, L to t = Lx + £ 7 , indicated by *. Several notable pulsars 
have been labeled. Blue squares: gamma-ray-selected pulsars. Red triangles: millisecond gamma-ray pulsars. Green 
circles: all other radio loud gamma-ray pulsars. Unfilled markers indicate pulsars for which only a DM-based distance 
estimate is available. Pulsars with two distance estimates have two markers connected with dashed error bars. 



sars. Among the other methods, the most commonly 
used estimates the pulsar distance from the radio dis- 
persion measure (DM) coupled to an electron density 
distri bution model, as the NE2001 Cordes and Lazid 
2002]. The distances measured with this method have 
to be taken carefully and a conservative error of 30% 
is assumed for the calculation of the pulsar luminosity. 
When distances from different methods disagree and 
no one is more convincing than the other, a distance 
range is assumed. 

In figure 2] the pulsar distribution on the Galactic 
plane clearly shows that the LAT is detecting a local 
population of gamma-ray pulsars, most of them not 
farther than 3 kpc (the inner circle in the figure). For 
the pulsars with an estimated distance range, only the 
closer edge of the range is drawn in the figure. 

In order to give a complete description of the 
gamma-ray pulsar population, the possible association 
with sources in other wavelengths has been investi- 
gated. At least 19 of the 46 LAT gamma-ray pulsars 
are positional associated with a PWN and/or SNR, 
and 9 of them are associated with TeV sources. Those 
pulsars with both TeV and PWN associations are typ- 
ically young, with ages less than 20 kyr. There are at 
least 3 more pulsars associated with TeV sources but 
not with PWNs or SNR. 

Among the LAT gamma-ray pulsars, 17 were 
unidentified EGRET sources. This result improves 



the discussions about the nature of the Galactic 
unidentified EGRET sources, as well as the discus- 
sion on the pulsars emission models. Population syn- 
thesis studies have shown that the expected ratio of 
radio-loud and radio-quiet pulsars strongly depends 
on the assumed emission model. The total Galac- 
tic birthrate of energetic pulsars evaluated from the 
LAT sample is of 1/50 yr, with gamma-selected ob- 
ject representing half or more. The small ratio of 
radio-selected to gamma-selected gamma-ray pulsars 
suggests that gamma-ray emission has an appreciably 
larger extent than the radio beams, such as expected 
in th e outer gap (OG) [Romani and Yadigaroglul 
1995 1 and slot-gap/two pole caust ic (SG/TPC) mod- 
els Muslimov and Harding! |2004| jDvks and Rudakl 
|2003| . 

Other constraints on pulsar models can come from 
the investigation of the relation between the pulsar 
luminosity (£ 7 ) and the rotational energy loss rate 
(E). In figure [5] is plotted this relation for the gamma 
ray pulsars detected by LAT. The dashed line indi- 



cates L 7 



E, while the dot-dashed line indicates 



oc 



E 1 / 2 , where 



L 7 = And /nGioo- 



(3) 



The flux correction factor /q Watters et al\\206^ 
is model-dependent and depends on the magnetic in- 
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clination and observer angles a and £. Both the 
outer gap and slot gap models predict /□ ~ 1, in 
contrast to earlier use of jfa = l/4n w 0.08 [in e.g. 
Thompson et aZ.|[l 994l. or / n = 0.5 for MSPs [in e.g. 



Fierro et all 1995^ For simplicity, we use fn = 1 



throughout the paper, which presumably induces an 
artificial spread in the quoted L 7 values. However, 
it is the quadratic distance dependence for L 7 that 
dominates the uncertainty in L 7 in nearly all cases. 

Gamma-rays dominate the total power £ to t radi- 
ated by most known high-energy pulsars, that is, 
L tot « Lj. The Crab is a notable exception, with 
X-ray luminosity Lx ~ 10L 7 . In Figure [5] we plot 
both L 7 and L x +L^. L x for E < 100 MeV is taken 



from Figure 9 of lKuiper et all 200 1| 

Many models for the gamma-ray emitting gap pre- 
dict that the gamma-ray luminosity of the pulsars 
is proportional to the root square of the spin-down 
power (L-y oc E). This imply that the gamma- 
ray efficiency increases with decreasing , down to 
~ 10 34 — 10 35 erg s _1 where the gap saturates at 
large efficiency. However the luminosity evolution 
of the LAT pulsars sample seems more complicated. 
At high E (> 10 36 - 5 erg s~ 4 ) the data seems to be 
in rough agreement with the predicted trend. For 
10 35 < E < 10 36 - 5 erg s~ 4 L 1 seem nearly constant, 
but the uncertainties in this range are very large and 
could hide the true luminosity evolution. In the range 



10 < E < 10 erg s the gap saturation is ex- 

10 35 erg s _1 ) and 



10 34 erg s x ) and at least for 



pected to occur in both Slot Gap ( 
Outer Gap models ( 
the OG model f is expect to drop well below 1. How- 
ever, with the wide range of gamma-ray efficiencies 
in the plot it is not possible to discriminate between 
these model predictions. Finally for E < 10 34 erg s _1 
the sample is dominated by the MSPs and the trend 
seems more consistent with E. 



4. DISCUSSION 

The LAT is providing a new local population of 
gamma-ray pulsars, composed by tree families: the ra- 
dio loud normal pulsars, the gamma-selected pulsars, 
and the millisecond pulsars. Considering also that 
gamma-ray pulsations have been searched for pulsars 
spanning the entire P — P plane, it makes of the LAT 
population an unbiased one that provide a new win- 
dow on pulsar demographics and physics. 

Evidence is that a large fraction of the local ener- 
getic pulsars are GeV emitters. There is also a sig- 
nificant correlation with X-ray and TeV bright pul- 
sar wind nebulae. Conversely, the pulsar origin of 
a large fraction of the bright unidentified Galactic 
EGRET sources is now uncovered. Plausible pulsar 
counterparts are also found for several previously de- 
tected TeV sources. In this sense the "mystery" of the 
unidentified EGRET sources is largely solved. 



Finally, the light curve and spectral evidence sum- 
marized above suggests that these pulsars have high 
altitude emission zones whose fan-like beams scan over 
a large portion of the celestial sphere. This is also con- 
firmed by the estimated pulsars birth rate and by the 
ratio of radio-selected to gamma-selected gamma-ray 
pulsars. 



Acknowledgments 



The Fermi LAT Collaboration acknowledges the 
generous support of a number of agencies and insti- 
tutes that have supported the Fermi LAT Collabo- 
ration. These include the National Aeronautics and 
Space Administration and the Department of Energy 
in the United States, the Commissariat a l'Energie 
Atomique and the Centre National de la Recherche 
Scientifique / Institut National de Physique Nucleaire 
et de Physique des Particules in France, the Agenzia 
Spaziale Italiana and the Istituto Nazionale di Fisica 
Nuclcare in Italy, the Ministry of Education, Cul- 
ture, Sports, Science and Technology (MEXT), High 
Energy Accelerator Research Organization (KEK) 
and Japan Aerospace Exploration Agency (JAXA) in 
Japan, and the K. A. Wallenberg Foundation and the 
Swedish National Space Board in Sweden. 

Additional support for science analysis during the 
operations phase is gratefully acknowledged from the 
Istituto Nazionale di Astrofisica in Italy and the Cen- 
tre National d'Etudes Spatiales in France. 

The Parkes radio telescope is part of the Australia 
Telescope which is funded by the Commonwealth Gov- 
ernment for operation as a National Facility managed 
by CSIRO. The Green Bank Telescope is operated by 
the National Radio Astronomy Observatory, a facility 
of the National Science Foundation operated under co- 
operative agreement by Associated Universities, Inc. 
The Arecibo Observatory is part of the National As- 
tronomy and Ionosphere Center (NAIC), a national 
research center operated by Cornell University under 
a cooperative agreement with the National Science 
Foundation. The Nancay Radio Observatory is op- 
erated by the Paris Observatory, associated with the 
French Centre National de la Recherche Scientifique 
(CNRS). The Lovell Telescope is owned and operated 
by the University of Manchester as part of the Jo- 
drell Bank Centre for Astrophysics with support from 
the Science and Technology Facilities Council of the 
United Kingdom. The Westerbork Synthesis Radio 
Telescope is operated by Netherlands Foundation for 
Radio Astronomy, ASTRON. 



eConf C091122 



6 



2009 Fermi Symposium, Washington, B.C., Nov. 2-5 



References 

Abdo, A. A., M. Ackermann, M. Ajello, B. Anderson, 
W. B. Atwood, M. Axelsson, L. Baldini, J. Ballet, 
G. Barbiellini, M. G. Baring, D. Bastieri, B. M. 
Baughman, et al, 2009a, Science 325, 840. 

Abdo, A. A., M. Ackermann, M. Ajello, B. Anderson, 
W. B. Atwood, M. Axelsson, L. Baldini, J. Ballet, 
G. Barbiellini, M. G. Baring, D. Bastieri, B. M. 
Baughman, et al, 2009b, ApJ Submitted (PSR 
J1907+0602 MGRO source). 

Abdo, A. A., M. Ackermann, M. Ajello, W. B. At- 
wood, M. Axelsson, L. Baldini, J. Ballet, G. Barbi- 
ellini, M. G. Baring, D. Bastieri, B. M. Baughman, 
K. Bcchtol, et al, 2009c, Science 325, 848. 

Abdo, A. A., M. Ackermann, W. B. Atwood, L. Bal- 
dini, J. Ballet, G. Barbiellini, M. G. Baring, 
D. Bastieri, B. M. Baughman, K. Bechtol, R. Bel- 
lazzini, B. Berenji, et al, 2008, Science 322, 1218. 

Abdo, A. A., et al., for the Fermi LAT collaboration, 
2009, ArXiv e-prints 0910.1608. 

Camilo, F., P. S. Ray, S. M. Ransom, M. Burgay, T. J. 
Johnson, M. Kerr, E. V. Gotthelf, J. P. Halpern, 
J. Reynolds, R. W. Romani, P. Demorest, S. John- 
ston, et al, 2009, ApJ 705, 1. 

Cordes, J. M., and T. J. W. Lazio, 2002, ArXiv 
e-prints (arXiv:astro-ph/0207156), arXiv:astro- 
ph/0207156. 

Dyks, J., and B. Rudak, 2003, ApJ 598, 1201. 



Fierro, J. M., Z. Arzoumanian, M. Bailcs, J. F. 

Bell, D. L. Bcrtsch, K. T. S. Brazier, J. Chiang, 

N. D'Amico, B. L. Dingus, J. A. Esposito, C. E. 

Fichtcl, R. C. Hartman, et al, 1995, ApJ 447, 807. 
Harding, A. K., J. V. Stern, J. Dyks, and M. Frack- 

owiak, 2008, ApJ 680, 1378. 
Hirotani, K., 2008, ApJL 688, L25. 
Kuipcr, L., W. Hermsen, G. Cusumano, R. Dichl, 

V. Schonfclder, A. Strong, K. Bennett, and M. L. 

McConnell, 2001, A&A 378, 918. 
Muslimov, A. G., and A. K. Harding, 2004, ApJ 606, 

1143. 

Ray, P. S., et al, 2009, ApJ In prep (Precise Timing 

of Fermi Gamma- Ray Pulsars). 
Romani, R. W., 1996, ApJ 470, 469. 
Romani, R. W., and I. -A. Yadigaroglu, 1995, ApJ 

438, 314. 

Takata, J., and H.-K. Chang, 2007, ApJ 670, 677. 
Thompson, D. J., Z. Arzoumanian, D. L. Bcrtsch, 

K. T. S. Brazier, J. Chiang, N. D'Amico, B. L. 

Dingus, J. A. Esposito, J. M. Fierro, C. E. Fich- 

tel, R. C. Hartman, S. D. Hunter, et al, 1994, ApJ 

436, 229. 

Watters, K. P., R. W. Romani, P. Weltevrede, and 

S. Johnston, 2009, ApJ 695, 1289. 
Zhang, L., K. S. Cheng, Z. J. Jiang, and P. Leung, 

2004, ApJ 604, 317. 



eConf C091122 



